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\  ABSTRACT  ft  ** 

^  This  note  obtains  the  theoretical  autocorrelation  function  of  an^ARMli) 
model  with  multiplicative  seasonality.  It  is  shown  that  this  function  can  be 
interpretated  as  the  result  of  the  interaction  between  the  seasonal  and 
regular  autocorrelation  patterns  of  the  ARMA  model.  The  use  of  this  result 
makes  easier  the  identification  of  the  structure  of  the  model,  is  helpful  in 
choosing  between  a  multiplicative  or  additive  seasonal  component  and  leads  to 
a  better  understanding  of  the  properties  of  the  estimated  autocorrelation 
function  of  scalar  ARMA  processes. 
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SIGNIFICANCE  AND  EXPLANATION 


In  the  process  of  building  an  ARIMA  model  for  a  time  series,  an  initial 
model  should  be  identified  analyzing  the  patterns  of  the  estimated 
autocorrelation  function  and  partial  autocorrelation  function.  Comparing 
these  observed  functions  to  the  theoretical  ones  associated  with  different 
ARMA  models  an  initial  model  can  be  entertained. 

This  identification  stage  was  sometimes  very  difficult  for  seasonal 
models,  because  the  theoretical  structure  of  the  autocorrelation  function  was 
not  completely  known.  This  function  is  obtained  in  this  paper.  The  result  is 
not  only  interesting  from  a  theoretical  point  of  view  but  has  important 
practical  indications  as  shown  in  an  example. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


THE  AUTOCORRELATION  FUNCTION  OF  SEASONAL  ARMA  MODELS 


Daniel  Pena* 

1.  INTRODUCTION 

It  has  often  been  stressed  that  the  most  difficult  task  in  obtaining  an  ARIMA  model 
for  a  given  time  series  is  the  identification  of  the  order  of  the  process.  Box  and 
Jenkins  (1970)  developed  broadly  the  theoretical  properties  of  the  autocorrelation 
function  (acf )  and  partial  autocorrelation  function  (paf )  for  processes  without  seasonal 
structure  and  outlined  the  pattern  of  the  acf  in  some  special  seasonal  processes  when  the 
regular  part  is  moving  average  of  order  one  or  two.  Cleveland  (1972)  proposed  the  inverse 
autocorrelation  function  as  an  alternative  to  the  paf .  Hamilton  and  Watts  (1978)  obtained 
the  exact  paf  for  the  simplest  ARIMA  models  and  were  able  to  show  the  general  pattern  of 
this  function,  as  a  simple  composite  of  the  autocorrelation  and  partial  autocorrelation 
coefficients  of  the  regular  component.  Finally,  Cleveland  and  Tiao  (1979)  have  proposed  a 
broad  class  of  seasonal  models  in  which  the  dependence  structure  among  the  observation  is 
not  invariant  to  shifts  in  time,  as  assumed  in  the  standard  seasonal  ARIMA  representation. 

One  important  difficulty  in  the  identification  of  seasonal  ARIMA  processes  is  that 
the  ignorance  of  the  exact  properties  of  the  acf  made  this  stage,  according  to  Hamilton 
and  Watts  (1978),  "a  perplexing  task  because  the  correlation  function  is  susceptible  to 
confusing  distortion  in  the  case  of  seasonal  time  series”.  We  will  show  in  this  note  that 
the  acf  of  a  seasonal  process  can  be  interpreted  readily  as  the  result  of  the  interaction 
between  the  seasonal  and  regular  components  of  the  model. 


^ ^ * 

*Daniel  Pena  is  Professor  of  Statistics  at  the  Escuela  Tecnica  Superior  de  Ingenieros 
Industriales,  University  of  Madrid.  The  author  is  indebted  to  Arthur  B.  Treadway  for  many 
helpful  comsients  on  an  earlier  draft  of  this  paper. 


Sponsored  partly  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041  and  the 
United  States-Spanish  Joint  Committee  for  Educational  and  Cultural  Affairs. 
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B: 


Let  ua  writ*  •  stationary  seasonal  Ann  eodel  aai 

Wfc  -  *(B>  ^(8*)  afc 

where  ♦  and  #a  ara  possible  infinite  operators.  Let  X(B)  and  (^(B*)  be  the 
autocovariance  generating  functions  for  the  regular  and  seaaonal  parta  of  Then,  it 

can  be  shown  (Box  and  Jenkins  (1970)),  that  the  autocovariance  generating  function  of  the 
process,  Y (B) ,  isi 

m 

Y(B)  -  0*  1(B)  ra(Ba)  -  l  TjB1  (2.1) 

2  ~“ 
where  is  the  variance  of  the  white  noise  process,  at,  and 

X(B)  -  *<B)  *(B_1)  -  l  X^1 
rs(B*)  -  *6(b“)  *#(b“*)  -  l  rgiB*i 


The  autocorrelation  generating  function  of  Wfc  is: 


a(B)  -  y-1  T(B)  -  o l  y;,X0r0P(B)Pa(B') 


(2.2) 


where  P(B)  ”  X  1X(B)  and  p  (B* )  -  r  'r  (B* )  are  the  autocorrelation  generating  functions 
o  I  ua 

of  the  regular  and  seasonal  part* 

Let  us  call  rA,  R1  the  theoretical  autocorrelation  coefficients  of  the  regular  and 
seasonal  part  of  the  nodel.  Then,  w*  can  write 

a(B)  -  Kp(B)pb<b“)  (2.3) 


where 


-1  2 .  _ 
Y_  a  x.r. 
0  a  0  0 


(1  +  2  l 


i-1 


r.iRsi) 


-1 


(2.4) 
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KB)  -  l  PhBh  -  K(  l  rjBjH  l  R^B**)  -  K  l  l  r^B1 


(2.5) 


•o  that  fay  using  the  value  of  K  and  equating  powers  of  B, 


r<  +  E  R.i  {r.i+j  +  r.i-j) 


i-t 


(2.6) 


1  +  2  T  r  .  R  . 
i-1  8i  8i 


Equation  (2.6)  displays  in  a  readily  understandable  way  the  structure  of  the  acf  for 
multiplicative  seasonal  models.  As  far  as  the  interpretation  of  (2.6)  is  concerned,  there 
are  two  different  cases.  The  first  occurs  when  the  regular  structure  has  autocorrelation 
coefficients,  rA  ,  that  are  nearly  zero  for  i  >  s/2.  In  this  case  the  denominator  of 
(2.6)  is  approximately  unity  and: 


rj  +  j,  R.i(r.i-J  *  W 


(2.7) 


so  that  we  will  observe:  (a)  in  low  order  lags  (i  <  s/2),  the  exact  regular  pattern;  (b) 
in  seasonal  lags,  the  exact  seasonal  structures;  (c)  in  lags  that  are  near  multiples  of 
the  seasonal  periods  (i  »  is  t  h;  h<s/2)  the  reproduction  of  the  regular  structure, 
symetrically  at  both  sides  of  the  seasonal  periods. 

When  the  regular  autocorrelation  coefficients  do  not  vanish  even  approximately  for 
i  >  s/2,  distortion  can  be  expected  in  the  above  pattern.  The  problem  will  be  especially 
acute  if  the  seasonal  period  is  low,  for  example  for  quarterly  data,  and  when  the 
autoregressive  operator  has  one  root  near  the  unit  circle. 

3.  M  EXAMPLE 

In  order  to  illustrate  the  kind  of  information  that  can  be  obtained  by  the  use  of  the 
properties  of  the  interactions  between  the  regular  and  seasonal  structure,  we  will 
consider  the  ozone  data  that  has  been  widely  modeled  by  different  procedures.  Box  and 
Tiao  (1975)  fitted  a  (0,0,1)  *  (0,1,1)12  ARIMA  model  to  this  series.  Abraham  and  Box 
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(1978)  have  shown  how  this  sodel  could  be  improved  through  s  deterministic  seasonal 
modelling.  Cleveland  and  Tiao  (1979)  have  presented  another  useful  approach  for  these 
data. 

Table  1  shows  the  estimated  acf  and  paf  for  this  series  seasonally  differenced  for 
the  period  November ,  1955,  throu^i  November,  1969.  The  model  fitted  by  Box  and  Tiao 
(1975)  was: 

(1-B12)Z  -  ( 1+. 14B) ( 1-.89B12 )a  (3.1) 

(.08)  (.02) 

o2  -  0.974  Q(37)  -  36.9. 

Although  the  Ljung-Box  statistic  Q  does  not  reject  the  model,  the  acf  of  the 
residuals  shows  significant  values  at  lags  2,  22  and  24. 

Using  the  theoretical  properties  of  the  acf.  Table  1  strongly  suggests  that  the 
regular  part  is  autoregressive.  This  fact  is  clear  from  the  specific  structure  of  signs 
of  both  the  acf  and  pact  at  both  sides  of  seasonal  lags.  As  far  as  the  seasonal  structure 
is  concerned,  there  is  some  evidence  of  AR  structure  because:  (a)  the  values  of  the  acf 
at  lag  24  is  not  only  significant  but  the  pattern  of  signs  around  lag  24  suggests 
interactions  and  (b)  the  value  of  the  acf  at  lag  36  is  almost  significant  and  at  both 
sides  of  lag  36  we  find  the  expected  pattern  of  signs  for  a  negative  coefficient. 

The  simplest  hypothesis  for  the  seasonal  part  is  then  an  ARMA  (1,1)  with  negative 
autoregresaive  parameter.  The  two  first  coefficients  of  the  paf  suggest  an  AR(2)  for  the 
regular  component.  The  estimated  model  1st 

(1+.13B12)(1-.15B-.15B2)A1  Z  -  (1-.87B12)a  (3.2) 

(.08)  ( .08) ( .08) 

o2  -  0.907  Q(34)  -  23.8. 

The  F  statistic  to  test  the  reduction  of  variance  in  model  (3.2)  versus  model  (3.1) 
is  5.76  that  is  highly  significant  with  a  "  0.005.  Furthermore  the  acf  and  pacf  of  the 
residuals  do  not  produce  any  doubts  about  the  adequacy  of  model  (3.2). 
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